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ABSTRACT 

The rate of internally-driven evolution of galaxy discs is strongly affected by the life- 
times of the spiral patterns they support. Evolution is much faster if the spiral patterns 
are recurrent short-lived transients rather than long-lived, quasi-steady features. As 
rival theories are still advocated based on these two distinct hypotheses, I review the 
evidence that bears on the question of the lifetimes of spiral patterns in galaxies. 
Observational evidence from external galaxies is frustratingly inconclusive, but the 
velocity distribution in the solar neighbourhood is more consistent with the transient 
picture. I present simulations of galaxy models that have been proposed to support 
quasi-steady, two-arm spiral modes that in fact evolve quickly due to multi-arm insta- 
bilities. I also show that all simulations to date manifest short-lived patterns, despite 
claims to the contrary. Thus the transient hypothesis is favoured by both numerical 
results and the velocity distribution in the solar neighbourhood. 
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1 INTRODUCTION 

After bars, spirals are the most prominent features of galaxy 
discs. The question of whether the spirals are short- or long- 
lived features is of major importance to the internal evolu- 
tion of disc galaxies. In this paper, I evaluate the evidence 
that spirals are transient features that change within a few 
dynamical times, but are superseded by fresh patterns in a 
recurrent manner. 

We still lack a complete, and widely accepted, the- 
ory for the origin of spiral patterns in galaxie s. There is 
good evidence (e.g. iKormendv fc Normanfll979l l that many 
prominent spiral patterns are found in barred galaxies, or 
are driven by tida l interactions, or perh aps even DM halo 
substructure (e.g. iDubinski et adfeoosl ). There seems lit- 
tle d oubt that a tidally-driven s piral pattern evolve s rapidly 
(e.g. ISalo fc Laurikain"enlll993l ; fDobbs et aZ.ll2010h. but re- 



peated excitations may be possible (|Bvrd fc HowardllT993) . 
If bars are long-lived features, as is generally believed (fo r 
a dissenting view see iBournaud. Combes fc Semelml f20u5). 
spirals driven by b ars could also be long-lived p atterns. How- 
ever, simulations ()Sellwood fc Sparkd Il988h suggest that 
spirals and bars could have differing pattern speeds, which 
seemed to be supported by observational data (Scigar, Chor- 
ney & James 120031; iButa et al.ll2009l), although the latter 
group recently changed their minds ( Salo et w feoiO)! 

Although many spirals in galaxies could be driven re- 
sponses, the ubiquity of the spiral phenomenon suggests oth- 
ers are likely to be self-excited features of discs. In particular, 
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all driving agents can be excluded from A^-body simulations 
of isolated stellar discs, which continue to manifest spiral 
patterns that must be self-excited. A satisfactory theory to 
account for self-excited spirals has yet to emerge, despite 
decades of effort (see revi e ws by Toomrdl 19771 ; lAthanassoulal 
11984 iBertin fc Linlll996l; ISeHwoodll2010al) . While all agree 
spirals are gravitationally-driven variations in the surface 
density of the old stellar disc, there is no consensus even on 
the expected lifetimes of the patterns. 

Here I first explain why the duration of spiral features 
is an important issue, and then elaborate on the brief dis- 
cussion of the lifetimes of spirals presented bv Binney & 
Tremaine l2008l. hereafter BT08. p. 526). I evaluate four dif- 
ferent types of evidence that bear on the question whether 
spirals are short- or long-lived patterns. 



2 INTERNAL EVOLUTION OF GALAXY 
DISCS 

It is not obvious that the rate of evolution of galaxy discs is 
affected by the lifetimes of the patterns. It may seem that 
if the time-averaged amplitude and pitch angle of transient 
spirals does not differ much from those of a steady long-lived 
pattern, the rate of change in the distribution of angular 
momentum and other quantities would be similar. 

However, the rate of angular momentum transport by 
spiral waves is given by the rate at which wave action 
(Lyndcn-Bcll & Kalnais 1972) is transported at the group 
velocity jToomrdll969l) which can differ from the gravity 
torque by the advective Reynolds stress (BT08, Appendix 
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J), also known as "lorry transport." Thus angular momen- 
tum transport is substantially reduced for long-lived pat- 
terns that invoke fe ed-back via the long-wave branch of the 
dispersion relation (Markl [l977l ). where tie advective term 
actually reverses the sign of angular momentum transport 
l|Sellwoodll2010al l. Transient spirals do noil involve the long- 
wave branch, and therefore redistribute angular momentum 

more effe ctively. 

Also, iLvnden-Bell fc Kalnais! l| 19721 ) showed that stars 
are scattered by a slowly changing potential perturbation 
only near resonances. More precisely, a sp ral potential that 
grows and decays adiabatically, i. e. on a time-scale long com- 
pared with the orbital and epicyclic periods, will not cause 
a lasting change to a star's orbit. Lasting changes do arise 
where wave-particle interactions become important near the 
resonances. Stars experience secular changes through "surf- 
ing" on the potential variations at corotation or through 
a periodic forcing close to their epicyclic frequency at the 
Lindblad resonances. 

Since exchanges between the stars arid the wave take 
place at resonances, large changes induced by a long-lived 
pattern would be confined to a few narrow resonances, the 
most important of which, the inner Lindblad resonance, 
must be "sh i elded" to avoid fierc e damping of the spiral 
|Markl 1 1974 iBertin fe Lull 19961) . The frequency widths 
of resonances of a steady wave depend only on the am- 
plitude of the perturbation. But resonances are broader 
for time-dependent waves, and multiple, short-lived distur- 
bances having a range of pattern speeds c lange the angular 
momenta of stars over large parts of the disc. 

Any mechanism to excite spirals rel es on the extrac- 
tion of angular momentum from stars in 1 he inner disc and 
its dep osition on stars in the outer disc (Xvnden-Bell & 
Kalnais Tl972F Resonant stars of a finite-amplitude wave will 
change their orbits as a result of these ai gular momentum 
exchanges and a spiral wave must decay once the resonances, 
which are the sources a nd sinks that excite it, have been 
exhausted. iBertinl (|l983h found order unity angular momen- 
tum changes would take a few Hubble times for a spiral of 
20% overdensity with a pitch-angle ~ 16 ' . But resonances 
will be depopulated on a much shorter tine-scale for large- 
amplitude, open spirals, which therefore < annot last long if 
the pattern speed is constant. Note that the lifetime could 
be lengthened by a changing pattern sp ;ed, which would 
allow the resonances to sweep over the dii c. 

Thus the question of the lifetimes of spirals is of ma- 
jor importance to our understanding of iisc evolution. If, 
on the one hand, the recurrent transient sicture is correct, 
spirals are the most important agents driying the evolution 
of galaxy discs: they transport angular momentum, scat- 
ter stars into non-circular orbits, and cause radial mixing 

erns, on the other 
hand, which must have moderate amplitude and be tightly 



wrapped in order to persist, would have a 
cant effect. 



3 OBSERVATIONAL EVIDENCE 



Near - IR photometric images ([Schweizerl 



much less signifi- 



976; 



Block et 

1994 iGnedin. Goodman fc Freilll995l; iRixl fc Zaritskvll 19951: 
Seigar fc James!ll998l ; lEskridge et al.ll2OO2l ; lGrosb0l. Patsis fc 



Pompe|2004) reveal that spiral patterns represent large am- 
plitude variations in the surface b rightness of the old stellar 
disc. lzi jetti. Chariot fc Rixl (|2009h attempt to derive a more 
faithful mass map from their H-band images by applying a 
pixel bj pixel correction based on population synthesis mod- 
els for the broad-band colours. All this work has confirmed 
that spiral patterns constitute large-amplitude density vari- 
ations i i the total stellar density in the disc, but says nothing 
about the lifetimes of the features. 

Ve' ocity maps of galaxies that can resolve inter-arm 
variatio as reveal that spiral patterns are associ ated with or- 
ganized non-axisymmetric streaming flows (e.g. IVisser|[l978l : 
ICanziai .. Allen & : Tilanuslll993t IShettv et qlteoOTT l . The am- 
plitude of the non-axisymmetric potenti al variations can be 
estimated by modeling such data (e.g. iKranz. Slvz fc Rixl 
2003), confirming that spirals are associated with significant 
potential variations. But unfortunately, these measurements 
are aga n quite insensitive to the lifetimes of the spirals; the 
flow pa ;tern is essentially the gas response to the instanta- 
neous p otential variation, and its qualitative features are in- 
depend ;nt of the mechanism that created it (BT08, p. 526). 
Since re sonances are broader for short-lived waves, the quan- 
titative gas response may depend weakly on the duration of 
the wa r e, but any attempt to use this as a discriminant 
betweei . theories would require an accurate independent es- 

timate )f the spiral amplitude. 

Gn adin. Goodman fc Frei (Il995l) and Fovle. Rix & Zi- 
betti l|2hlfF estimate the gravity torque from the spiral pat- 
tern. If this measurement determines the rate of angular 
momen mm transport, then one can use it to estimate the 
time sc lie for significant changes to the initial angular mo 
mentur i distribution. If the spiral is long-lived, the grav- 
ity tore ue must be corrected for the advective term (32J), 
which can substantially lengthen the time-scale. Advection 
is far le is important for transient spirals. While it is of great 
interest that the angular momentum evolution time scale 
could be less than a Hubble time were spirals to be tran- 
sient, tl is measurement again does not discriminate between 
short- c nd long-lived spiral patterns. 

Thp "density wave" theory of spiral modes (33]3) re- 



quires t 
Toomrc 
range 1 
a larger 



lenging 
which i; 



he outer disc to be dynamically cool, in particular 
s axisymmetric stability parameter should be in the 
< Q < 1.2. In principle, a direct measurement of 
value of Q would cause difficulties for this theory. 
A num )er of attempts have been made to est imate Q in 
real die zs ( e.g. lKormendvlll984 iBottemal Il993l; Herrmann 
Ciardullo 120091) . but uncertainties are large because the 
measur ;ment of the low velocity dispersion in a disc is chal- 
and the Q value depends on the disc surface density, 
i not tightly constrained. It should also be noted that 
the raw] Q-value estimated from the stars only may not be 
the relevant quantity in realistic discs, where the gas com- 
ponent and perhaps also the young stellar population can 
have a l isprop ortionate effect on the overall responsiveness 
dRafiko l[200lh- 

Doabs & Pringle (|2010h suggest that the downstream 
age disi ribution of star clusters can be used to determine 
the origin of a spiral pattern. However, their method re- 
test only of whether the gas that formed the star 
streamed through the spiral arm at a speed that dif- 



ally is 
clusters 



whethe 



fers significantly from the local circular speed - i.e. merely 



a coherent density wave is present - and does not 
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discriminate between mechanisms for its origin or test for 

its lifetime. 

Meidt et al. ll200^) generalize the method of Trcmaine 
& Weinberg (|1984l ) to discs that may have multip le pat- 
tern speeds. They (|Meidt. Rand fc Merrifieldl l2009h apply 
it to CO data on a number of galaxies, finding evidence for 
multiple patterns in several cases. Other methods to esti- 
mate the pattern speed of spiral arms often make the (im- 
plicit) assumption that there is a single pattern; for example , 
the kinematics of stars formed in spiral arms (|Bashlll979l) . 
the change in radial flow direction across corotation (e.g. 



lElmegreen. Wilcots fc Pisancj 19981). or th e change in the 



pattern of resid ual velocities Icanzianll 19931) , or brute force 
modelling (e.g. iRautiainen et al. 20081) . etc. Once again, 
however, the results tell us little about the lifetimes of the 
patterns. 

In summary, I concur with BT08 that it is frustratingly 
difficult to determine the lifetimes of spiral patterns, or the 
mechanism that gave rise to them, from observations of ex- 
ternal galaxies. 

However, evidence from the Milky Way seems to favour 
the transient spiral picture. The in-plane velocity compo- 
nents of stars in the s olar neighbourhood as revealed by the 
HIPPARCOS mission (|Dehnenlll99Sl; iNordstrom et qZ.ll20o4) 
do not have the smooth double Gaussian anticipated by K. 
Schwarzschild (BT08, p. 321), but instead are characterized 
by a number of separate str eams, with essentially no un - 
derlying smooth component (|Bovv. Hogg fc Roweisl f2009). 
The features are too substantial to have simply arisen from 
groups of sta rs that were born with similar kinematics (e.g. 
Egge nl 1 19961). as con firmed in detailed studies (Famaey et 
al. l2007t iBensbv et aZ.ll2007l ; iBovv fc Hogell20ia ). and it is 
clear that the entire DF has been sculptured by dynamical 
processes. 

Individual 
tributed to resonances with the bar 



features in this distribution have been at- 

~oi1ll998l; 



(jRaboud et 

lDehnenll2000l;|Fuxll200ll) or with a spiral pattern (Quillen & 
Minchev Bo05l ; ISellwoodll2010bl), while other models inc lude 
both bars and spirals (lQuillenll2003l; IChakrabartvll2007l; An- 
toja et al. 120091 ). The overall appearance of the velocity dis- 
tribution was successfully reproduced by De Simone, Wu 
& Tremaine (|2004l ) who invoked a succession of short-lived 
spiral transients. 



4 SPIRAL STRUCTURE THEORY 

Three distinct mechanisms have been proposed to account 
for self-excited spiral patterns in disc galaxies. Bertin & 
Lin ( 1996) propose that spiral features are manifestations of 
quasi-steady global modes of the underlying disc. Goldreich 
& Lynden-Bcll (19651) and IToomrd (|l990h argue that, aside 
from tidal- and bar-driven cases, most spirals are short-lived 
transient features that are essentially collec tive responses to 
internal density fluctuations within the disc. Scllwood (1 20001 ) 
also supposes that spirals are short-lived, but suggests they 
result from a recurrent cycle of vigorous l arge-scale modes . 
Here I give brief summaries of each in turn; ISellwoodl (|2010aD 
gives a more detailed review. 



4.1 Swing Amplified Noise 

Swing amplificatio n was first described by Go ld reich fc 
Lynde n-Bell (|l965l ) and I Julian fc Toomrd l|l966l) . IToomrd 
(|l98li) highlighted its importance, which he illustrated with 
the now classic figure named "dust-to-ashes" reproduced 
in BT08 (their Fig. 6.19). It shows the dramatic tran- 
sient trailing spiral that results from a small input leading 
disturbance. lJulian fc Toomrel (| 19661) showed that swing- 
amplification causes the collective response of a disc to a 
co-orb iting ma ss clump to be a substantial spiral "w ake." 

Goldreich & Lynden-Bell i| 19651) and IToomrd (|l990l) 
suggested that a large part of the spiral activity observed in 
disc galaxies is the collective response of the disc to clumps 
in the density distribution. As a spiral wake is the collective 
response of a disc to an individual co-orbiting perturber, 
multiple perturbers will create multiple responses that all 
orbit at different rates. The behaviour of this polarized disc 
reveals a changing pattern of trailing spirals, which can 
equivalcntlv be regarded as swing-amplified noise. Toomre 
& Kalnajs ( 199 ll ) show that amplified noise arising from the 
massive disc particles themselves can be understood in the 
shearing sheet, where the resulting spiral amplitudes are lin- 
early proportional to the input level of shot noise. 

This could be a mech anism for "a swirling hotch-potch 
of pieces of spiral arms" (|Goldreich fc Lvnden-Belilll965l ) 
in very gas rich discs, where a high rate of dissipat ion may 
be ab le to maintain the responsiveness of the disc (|Toomrd 
1990) while the dumpiness of the gas distribution may make 
the seed noise amplitude unusually high. However, the spiral 
structure should be chaotic, with little in the way of clear 
symmetry expected. Also, it seems likel y that spiral ampli- 
tudes (e.g. IZibetti. Chariot fc Rixi 2009) are too large to be 
accounted for by this mechanism in most galaxies. 



4.2 Recurrent Cycle of Groove Modes 

Sellwood & Kahn (|l99ll) showed that discs are destabilized 
by a deficiency of stars over a narrow range of angular mo- 
mentum; in a disc without random motion, such a deficiency 
would be a "groove." The groove itself is unstable, and the 
instability becomes a global mode through the vigorous sup- 
porting response of th e surrounding disc , prod ucing a large- 
scale spiral instability. ISellwood fc Kahnl (|l99lt ) showed that 
the linear instability, which is driven from corotation, also 
develo ps in discs with ra ndom motion, while Sellwood & 



Binnev 1[2002Tr showed that the amplitude of the mode was 
limited by the onset of h o rsesho e orbits. 

Sell wood fc Kahn Jl99ll), and earlier Lovelace & 
Hohlfeld (fl978h . found that almost any narrow feature in 
the angular momentum density is destabilizing. Thus the 
common starting assumption of spiral structure studies, that 
the underlying disc is featureless and smooth, may throw the 
spiral baby out with the b athwater. 

Sellwood & Lin (|l989l ) showed, in simulations of a low 
mass disc, that the particles were scattered at the Lindblad 
resonances as each coherent wave decayed, thereby leaving 
behind an altered distribution function that cr eated the con- 
ditions for a new instability. ISellwoodl (|2000l ) described in 
more detail that each spiral pattern seems to be a strongly 
unstable true mode of the disc in which it grows; it lasts 
for only a few rotations, changing the disc properties as it 
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decays so as to provoke new instabilities of the altered disc 
having different pattern speeds. The discovery of resonance 
scattering among the stars of the solar neighbourhood (Sell- 
wood |2010bJ) suggests that this mechanism may occur in 
real galaxies. I am actively engaged in addressing the many 
details of this picture that remain unclear. 

4.3 Long-lived Global Modes 

Simple models of disc galaxies support ma ny global linear in- 
stabilities (e.g. lKalnaislll978l ; |jalalill2007h . The bar-forming 
mode is generally the fastest growing, but has almost no spi- 
rality. These studies are therefore important to understand 
stability, but do not appear promising for spiral generation. 

The "de nsity wave" th e ory fo r spiral modes, described 
in detail by iBertin fc Linl fl996), invokes a more specific 
galaxy model with a cool outer disc and h ot inner disc. T he 
local stability parameter, Q = <jr/ 'oR }C rit l|Toomrelll964h . is 
postulated to be Q > 1 in the outer disc and to rise steeply 
to Q > 2 near the centre. Under these specific conditions, 
Bertin & Lin fi nd slowly e volving spiral modes that grow by 
their WASER l|MarkHl977l ) mechanism. They invoke shocks 
in the gas to limit the amplitude of the slowly gr owing mode, 
leadin g to a quasi-steady global spiral pattern. lLowe et al\ 
jl994) present a model of this kind to account for the spiral 
structure of M81. 

The principal objection to their picture is that it is likely 
that an outer disc with such a low Q will support other, 
more vigorous, collective responses that will quickly alter 
the background state by heating the outer disc, as I now 
show. 

4.4 A Direct Test 

Here I report simulations designed to te st directly whether 
a galaxy model of the kind invoked by (|Bertin et qZ.I IT98Sl . 
hereafter BLLT) can in fact survive to support the slowly- 
growing global mode they predict should dominate. 

The model I adopt here is from the Ai survey of BLLT, 
which is a minor variant of that p roposed by Fall & Efs- 
tathiou (|l980f v" The surface density of the disc is 

Z(R) = (l + A)^e R / R «f(2R/R d ), (1) 

with 

/(2/H (l-(l+%)(l-y) 4 (l-e-«/^/6) V<1 (; 
I- 1 otherwise. 

Here Rd and Md are respectively the scale length and mass, 
when A = 0, of the unmodified exponential disc. The f(y) 
factor creates a central dip in the surface density of approx- 
imate radius Rd/2. The rotation curve of their model has 
the form 

V(r) = Vb — fL— with x = 2R/R d , (3) 

that is independent of the disc contribution. Because mass 
removed by the central cutout is effectively replaced by rigid 
matter, the model can be thought of as including a small 
central bulge as well as a pseudo-isothermal halo. Here I 
adopt units such that G — Md ~ Rd = 1, and also Vb = 1; 
the orbit period at 7? = 2Rd is ~ 13 time units. 



1.2 




1 1 1 1 L 

2 4 

Radius 



Figure 1. Top panel shows the rotation curve of the model dis- 
cussed in this section in units given in the text. The solid line 
is the full circular speed, the dotted line shows the contribution 
from the modified exponential disc, and the dot-dashed line that 
from the halo. The lower panel shows the initial radial variation 
of Q. 



BLLT also specify the radial variation of Toomre's lo- 
cal axisymmetric stability parameter Q(R) = Qod(1 + 
l.5e 2R ^ Rd ); i.e. Q = 2.5Qod at the centre and decreases 
rapidly to Q = Qoo over most of the outer disc. Be- 
cause BLLT do not supply stellar dynamical distribution 
functions, I set the initial velocities to create a rotation- 
ally supported disc having the desired Q profile, using the 
Jeans equations in the epicyclic approximation to deter- 
mine the azimuthal dispersion and asymmetric drift (BT08, 
eq. 4.228). 

I have chosen the case with Qod = 1 and A = —0.35. 
The upper panel of Fig. [1] shows the rotation curve, which 
indicates that this model has a strongly sub-maximal disc 
l|Sacketj|l997f ). i.e. the disc contributes merely ~ 25% of 
the central attraction at R = 2R d . The lower panel shows 
the initial Q-profile. Because the disc is so light, and the 
disc surface density dips towards the centre, the initial ran- 
dom velocities are small enough that the epicycle approxi- 
mation is adequate and the Jeans equations lead to a quite 
respectable equilibrium model. 

BLLT predict that this model has the slowly-growing, 
tightly-wrapped, bi-symmetric spiral mode illustrated sec- 
ond from left in the bottom row of their Fig. 3, but they do 
not give the frequency of the mode. 

Here I report the results of two series of simulations 
with this disc model. In the first series, I restrict disturbance 
forces torn = 2 only, as do BLLT in their mode calculations. 
In the second, I allow sectoral harmonics 2 ^ m ^ 8 to con- 
tribute. In neither case do I include disturbance forces from 
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Figure 2. Contours of power as a function of radius and fre- 
quency computed from the relative bi-symmetric density varia- 
tions (E2/S0) on the rings of the simulation grid. Peaks indicate 
coherent features with pattern speed mQ p that are present for a 
significant fraction of the simulation. The solid curve shows the 
circular angular frequency mf2, and the dashed curves m!i ± k. 
The simulation was run for 500 time units, or 38.5 disc rotations 
at R = 2. 

m = & 1. Excluding m = 1 avoids unbalanced disturbance 
forces that could arise from a lop-sided mass distribution in 
a rigid halo; there are no external perturbations to excite 
such waves, which are also most unlikely to be self-excited 
in this low mass disc (see below). The unchanging axisym- 
metric central attraction maintains the rotation curve shown 
in Fig. [T] 

The simulations employ the grid code that has been 
shown to reproduce global modes predicted by linear the- 
ory for a nu mber of stellar dynamical models ( Sellwood fe 
Athanassoulall98&ISeUwoodll989l ; lEarn fc Sellwoodl [l995l; 
ISellwood fc Evansll200lf) . The particles move over a 2D po- 
lar grid having 65 x 96 mesh points and their motion is 
integrated with a time step of 0.05. Direct tests reveal that 
results are insensitive to these numerical parameters. The 
particles interact with forces derived from a Plummer soft- 
ening kernel with scale 0.05i?d- 

With m = 2 only disturbance forces, a simulation with 
N = 200 000 particles barely evolved. No non-axisymmetric 
features were visible in the first 40 disc rotations and the 
Q-profile remained unchanged. There is a hint in the power 
spectrum shown in Fig. [2] of a coherent wave in the inner 
disc rotating at the angular rate mQ p ~ 0.8, but it is barely 
stronger than the noise peaks along the mQ. curve. With- 
out the predicted eigenfrequency, I am unable to confirm 
whether it corresponds to the WASER mode predicted by 
BLLT. Nevertheless, it is clear that the model does not pos- 
sess any rapidly-growing bi-symmetric instabilities. 

The outcome is completely different when disturbance 
forces from higher sectoral harmonics are included, as shown 
in Fig. [3] Irrespective of the number of particles, or whether 
the particles are positioned carefully on rings (a quiet start) 
or placed at random, the model quickly develops multi-arm 
spiral patterns that heat the disc, as shown in Fig. [4] 

The reason for this different behaviour when higher 
sectoral harmonics are included is clearly related to swing- 
amplification, which is strongly dependent on the value of 
X = Rk crit /m (| Julian fc Toomrell 196fj ; pToomrell 198 ll BT08), 




Figure 3. The time evolution of the simulation in which distur- 
bance forces included sectoral harmonics m > 2. Only one particle 
in 400 is plotted in order to avoid saturating the images. Multi- 
arm patterns appear soon after the start. The grid edge (outer 
circle) has radius R = 6.5 and one rotation at .R = 2 is 13 time 
units. 




Radius 



Figure 4. The radial variation of Q at three times in a simulation 
with 20M particles and sectoral harmonics 2 ^ m ^ 8 all active. 
The initial Q-profile is that shown in Fig. [T] 

where fc cr it = k 2 /(2tvGT,), the smallest wav enumber of ax- 
isymmetric Jeans instabilities (jToomrdll964T ). In the present 
model, selected from BLLT, this parameter rises steadily 
with radius from X ~ 8/m at R = 1 to X ~ 20/m at 
R = 3, continuing to higher values farther out. Since the 
swing-amplifier is vigorous only when 1 < X < 2.5 in a flat 
rotation curve, as here, it is pretty much dead at all radii 
for m — 2 waves, but vigorous responses are expected for 
m ^ 4, as we observe. 

Fig. [5] shows the time evolution of Q at R = 3Rd in 
three simulations with differing numbers of particles showing 
that heating of the disc, which coincides with the occurrence 
of visible spiral patterns, is increasingly delayed as larger 
numbers of particles are employed. The principal effect of 
increasing N is to reduce the amplitude of initial density 
fluctuations. The swing-amplifier quickly polarizes the disc, 
creating trailing spiral responses having an amplitude pro - 
portional to the input noise signal (Too mre fc Kalnaiall991fl . 
The larger the particle number the smaller the initial am- 
plitude, and it is clear from Fig. [5] that the swing-amplified 
noise has too small an amplitude to cause heating at first 
when N — 20M. We do not expect global instabilities in 
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Figure 5. The time evolution of Q at R = in three sim- 

ulations with differing numbers of particles. Sectoral harmonics 
2 ^ m ^ 8 were active in all three cases. 

a smooth disc for m > 2 because small-amplitude distur- 
bances at most reasonable patt ern speeds will be damped 
at an inner Lindblad resonance l|Marklll974[ ). Thus the later 
increase in the spiral amplitude must be a non- linear effect, 
pe rhaps related to the re current instability cycle reported 
bv lSellwood fc Linl |l989l) . Whatever the cause, the similar 
heating rate and final Q value is consistent with spiral am- 
plitudes that are independent of iV at later times. The origin 
of multi-arm spiral waves will be followed up in future work. 

Fig. [5] shows that even with N = 20M, heating begins 
after just 15 disc rotations. This should be contrasted with 
the case where forces were restricted to m = 2 where no 
visible bi-symmetric features appeared in 40 rotations even 
with the ten times larger seed amplitude implied by 100 
times fewer particles. It should be further noted that real 
spiral galaxy discs have large star clusters and giant molec- 
ular clouds that imply much larger seed density fluctuations 
than arise from iV = 20M randomly-placed, equal-mass par- 
ticles. 

4.5 Discussion 

While these simulations explicitly test just one of the many 
models presented in BLLT, all the cases that they regard 
of most interest have strongly sub-maximal discs in order 
that swing amplification is ineffective for m = 2; this aspect 
is essential so that the model possess only slowly-growing, 
tightly-wrapped, bi-symmetric spiral modes. However, the 
vigour of the swing-amplifier for m > 2 makes it inevitable 
that every one of their galaxy models with a low-mass disk 
will be subject to stronger activity due to disturbance forces 
with m > 2, which will quickly heat the outer disc and 
destroy the conditions they require, as just demonstrated. 
Thus the "basic state" invoked by BLLT for the modes they 
favour to account for bi-symmetric spiral patterns could not 
survive in real discs that permit disturbances of all sectoral 
harmonics. 

As this conclusion depends largely on the results from 
simulations, one must worry whether they can be trusted. 
The principal source of concern is that the origin of the 
multi-arm patterns that heat the disc remains obscure. Sim- 
ulations of this kind over many ye ars (e.g. Sellwood & Carl- 
berg ll984lFuiii et aZ.l2010h have manifested recurrent multi- 
arm spiral patterns, and the behaviour has not changed as 



numerical quality has risen and the codes have passed many 
tests. However, the possibility that the behaviour could re- 
sult from some artefact in the simulations cannot be ex- 
cluded altogether until a satisfactory explanation is pro- 
vided. Note that to doubt the result on these grounds calls 
into question all simulations of isolated discs over the past 
40 years, as well as those that model the formation of disc 
galaxies. 

Other possible criticisms, such as the simulations could 
be unreliable because of gravity softening, particle noise, or 
the restriction to 2D, are more readily rebutted. First, the 
very same code has been shown to reproduce modes pre- 
dicted from linear stability analyses of other models, as cited 
above. Also, Plummer softening in simulations with particles 
confined to a plane provides a reasonable allowance for fi- 
nite disc thickness and anyway the same beh aviour persists 
in simulations of discs with finite thickness (|Roskar et al\ 
l2008at iFuiii et aZ.ll2010h . Further, since the only variation 
in the behaviour as N is increased 100-fold is an increas- 
ing delay due to a decreasing seed amplitude, with no other 
qualitative differences, an argument that simulations cannot 
be trusted because N is too small is somewhat threadbare. 

One might also worry that the stellar dynamical real- 
ization I have created differs from the model that BLLT 
analyzed in the hydrodynamic approximation - essentially 
I have replaced pressure in their hydrodynamic calculations 
with velocity spreads; a more direct test would require a 
prediction of a slowly-growing spiral in a stellar-dynamical 
model. However, it is hard to see why this minor difference 
should matter, especially in this case where the dispersion 
is a small fraction of the orbit speed almost everywhere; the 
exception is at the centre where Q is so large that the disc 
is (by design) dynamically inert. 

The simulations here have not, of course, made any al- 
lowance for the influence of the gas component, which can 
offset the heating effec ts from transient spiral patt erns to 
some extent. However, ISellwood fc Carlberd (| 19841 ) found 
that simulations that included a very crude form of cooling 
still settled to 1.5 < Q < 2. 



5 SIMULATIONS 

Short-lived recurrent spiral patterns have developed spon- 
taneously i n simulations of isolated galaxies, fr om the first 
studies by iMiller^ rcndcraast & Quir kl (119701) and Hock- 
ney & Brownrigg l|l974h . right through to modern simula- 
tions that include more sophisticated physical process es (e.g. 
iRoskar et al.ll2008al ; lAgertz, Tevssier fc Moore1l2010h. 

Claims of long-lived spiral waves (e.g. iThomasson et all 

1990) have mostly b een based on simulations o f shor t du- 
ration. For example, lElmegreen fc Thomassonl (|l993l) pre- 
sented a simulation that displayed spiral patterns for ~ 10 
rotations, but the existence of some underlying long-lived 
wave is unclear because the pattern changed from snapshot 
to snapshot. Other claims are equally doubtful, as I show 
next. 

5.1 Direct Tests 

AslDonner fc Thomassonl <|l994l, hereafter DT94) and lZhand 
jl99d hereafter Z96) have presented evidence for long-lived 
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Figure 6. The rotation curve of the initial model adopted by 
DT94 and Z96 in units given in the text. The solid line is the full 
circular speed, the dotted line is the speed due to the modified 
exponential disc, and the dashed and dot-dashed lines show the 
speed due the bulge and halo respectively. 



spirals in the same model, I have chosen to try to reproduce 
their results here. I first summarize the model they employed 
and then report my own analysis of the similar results I 
obtain when I reproduce their simulations. 

DT94 adopted the disc surface density distribution 
l|Rohlfs fc Kreitschmarm1ll98fi l 



E(R) 



2M d 
2,-kR 2 , 



-R/Rd 



-2R/ R d ~\ 



(4) 



Here R d is the scale length of the outer exponential disc 
and M d is the disc mass. DT94 and Z96 chose M d = 0.5M t , 
where Mt is the total mass of the model, and employed two 
additional mass components to represent a central bulge 
and a halo, both of which exert the central attraction in 
the mid-plane of a razor-thin simple exponential disc. The 
masses and scale lengths were respectively 0.1M t , 0.1R d for 
the bulge and QAMt, 0.5R d for the halo. The rotation curve 
of this model is shown in Fig. [6] which compares well with 
that shown in Fig. 1 of DT94. These authors set the initial 
velocities in the disc such that Q = 1 at all radii. 

Here I recreate this model, and compute its evolution 
using essentially the same 2D polar grid code, but with a 
larger number of particles. The disc has iV = 2M particles 
that move over a grid having 100 x 128 mesh points. As 
DT94 and Z96, I use a Plummer softening law with a length 
scale 0.15R d to compute forces between particles. However, I 
adopt a more physically motivated set of units for which G = 
M t — Rd = 1. For comparison with the previous results, it 
should be noted that R d — 10 in their units, and one rotation 
period at R = 2R d , which takes 2nR/V c = 18.5 of my time 
units (925 time steps), is 314 time steps in DT94 and 628 
time steps in Z96. 

Fig. [7]shows the evolution computed here, which should 
be compared with that shown in Fig. 2 of Z96. Since the ten 
times larger number of particles used here lowers the seed 
amplitude, a little more evolution is needed for the spiral to 
grow. To make the closest possible comparison, I therefore 
show snapshots that are spaced at the same time interval, 
but are shifted later by a little over one disc rotation from 
the start. The overall appearance is quite similar; a strong 
m = 2 spiral is developing by time 72 and is perhaps more 
persistent than that in Zhang's calculation, where the spiral 
has faded more by the last two times. 




Figure 7. Evolution of the model run to reproduce that reported 
by DT94 and Z96. The six snapshots, which include one particle 
in 40, should be compared with those shown in Fig. 2 of both 
papers. The circle is drawn at R = 8R d , although the grid extends 
to 12.7R d . 



Since DT94 and Z96 claim that the m — 2 features are 
a long-lived spiral, I examine them more closely here. Fig. [8] 
shows that in my simulation they appear to be the super- 
position of several waves having differing pattern speeds. 
This figure should be compared with Fig. 5 of DT94, which 
presents a similar analysis for their model. 

Although I employed 40 times the number of particles 
used by DT94, the power spectrum is still quite noisy. How- 
ever, the lowest panel has three or more horizontal ridges 
that are caused by coherent waves extending roughly from 
the inner Lindblad resonance to a little outside the outer 
Lindblad resonance in each case. The lowest frequency peak, 
which is the farthest out in the disc, has the largest rela- 
tive amplitude, which is simply a reflection of the fact that 
the disc surface density decreases outward. Forming separate 
power spectra on the first and second halves of the evolution 
(top two panels of Fig. [5} reveals that the separate patterns 
reach peak amplitude in sequential order, with the fastest 
rotator (mfi p = 0.42) developing and decaying first - there 
is no significant power at that frequency in the second half of 
the run. At least two waves co-exist at significant amplitude 
for most of the evolution. 

A least-squares fit to these data ( Sellwood & Athanas- 
soula 19861 ). as well as to an expansion of the particle dis- 
tribution in logarithmic spirals, finds at least four coher- 
ent waves with mil p ~ 0.42, 0.30, 0.25 & 0.12, in rough 
agreement with the locations of the peaks in the lowest 
panel of Fig. [5] These values need to be multiplied by 
18.5/314 « 0.06 to convert to the frequency unit used by 
DT94, and these authors plot pattern speed fi p and not 
mQp that I measure; thus the pattern speeds I observe are 
0.013, 0.009, 0.007 & 0.004 in the units shown in Fig. 5 of 
DT94, where three waves can also be discerned. The fastest 
and most slowly rotating waves in DT94 have similar fre- 
quencies to those I observe, but neither intermediate wave 
corresponds to a feature in their plots. Exact correspondence 
should not be expected as all these waves are seeded by par- 
ticle noise. 

Fig. [9] confirms that the amplitude of bi-symmetric dis- 
turbances rises in the inner disc (dotted line) earlier than 
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radius 

Figure 8. Contours of power as a function of radius and fre- 
quency computed from the relative bi-symmetric density varia- 
tions (S2/S0) on the rings of the simulation grid. Peaks indicate 
coherent features with pattern speed mfl p that last for a signifi- 
cant fraction of the simulation. The solid curve shows the circular 
angular frequency mf2, and the dashed curves rail ± k. The top 
panel is computed from the first half of the simulation, the mid- 
dle panel from the second half, while the third panel shows the 
spectrum from both halves combined. 



in the outer disc (dot-dash line). Furthermore, the poten- 
tial variations fluctuate in amplitude, especially in the inner 
disc, in a manner characteristic of beats between waves ro- 
tating at different pattern speeds, as to be expected from 
the power spectra shown in Fig. [8] Were the evolution dom- 
inated by a single m = 2 spiral mode, the temporal evolution 
in this Figure would be smooth and vary synchronously at 
each radius. 

Thus "the bisymmetric spiral" in this simulation is not 
a single long-lived pattern, but the sup erposition of t hree, or 
more, waves that each grow and decay. [Zhand (|l998T ) argues 
that this model has a tendency to form a bar, which I have 
also noticed, and she therefore presents another simulation 
in which the disc mass is reduced to 0.4Af t and the halo 
mass increased to 0.5M t in order to weaken the tendency 
to form a bar. As she again claims that this new model 
supports a long-lived spiral, I have also tried to reproduce 
this slightly different simulation. As above, I find that the 
strong, and visually impressive, bi-symmetric disturbance is 
the superposition of more than two patterns that each do not 
last long. In this case the first wave to appear has frequency 
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Figure 9. The time variation of the relative amplitude of the 
bi-symmetric potential at three different radii. 



2f2 p « 0.29, but later the two dominant waves have 2Q P w 
0.36 & 0.19, with also some significant amplitude at m — 3. 



6 INDIRECT EVIDENCE 

Here I summarize three further arguments that bear indi- 
rectly on the question of spiral lifetimes. All three are related 
to the effect of spiral activity on the orbits of stars in the 
disc. 



6.1 Importance of Gas 

Transient spiral activity scatters disc stars away from circu- 
lar o r bits dBarbanis fe Woltierl Il967l ; ICarlberg fe Sellwood! 
1 19851: iBinnev fe Lacevll 19881) . As random motion of the disc 
stars rises, the collective behaviour that produces spiral pat- 
terns is weakened. Thus, in the absence of cooling, spiral 
activity is self-limiting - it heats the disc to a certain level 
at which the disc can no longer s upport spiral patterns. 

Sellwood & Carlberg (| 1984f) showed this explicitly in 
their simulations, and went on to show that a moderate 
amount of cooling in the form of fresh stars added to the 
disc on circular orbits enabled spiral activity to continue "in- 
definitely. "Subsec|uent - \vo^ & 
Freedman ll985l; lToomrelll990l ; I Roska r et al. I l2008al ), as well 
as modern galaxy formation simulations (Agertz, Teyssier 
& Moore 2010] and many others), has confirmed this re- 
sult in simulations with ever greater realism. Physically, the 
random motions gained by gas clouds are dissipated in colli- 
sions so that they keep moving on near circular orbits. Stars 
that form from the clouds therefore have similar motions, 
and a continuous supply of fresh stars on near-circular or- 
bits maintains a responsive stellar distribution that allows 
spiral activity to continue. Quantitatively, fresh stars added 
to a disc at the rate of a few every year is sufficient for this 
purpose. 

Thus the transient spiral picture offers a natural expla- 
nation for the absence of spiral patterns in SO disc galaxies 
that have little or no gas and the long-noted correspondence 
that disc galaxies with significant gas components also mani- 
fest spiral patterns. Note that this appealing argument again 
do es not un iquely favour short-lived waves, since Bertin & 
Lin (1 19961) also invoke a gas component both to limit the 
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amplitude of their mildly growing modes and to maintain 
the dynamically cool outer disc. 



6.2 Scattering of Stars 

It has been clear for some time that the velocity dispersion 
of disc stars in the solar neighbourhood rises with age (Wic- 
len ll977 l; lNordstrom et al$2004) and also, for main sequence 
stars, with colour ijAumer fc Binnevil2009f ). which is a surro- 
gate for mean age. T he highest ve l ocities cannot be produced 
by cloud scattering |Lacevlll99ll ; lHanninen fc Flvnnll2002l ) 
and some other accelerating agent, such as transient spirals, 
seems to be required. Long-lived spirals, which do not heat 
the disc nearly as effectively, could not achieve the requisite 
high velocities. 



6.3 Radial Mixing 

Studies of the metallicities and ages of nearby stars fEd- 



vardsson et al.\ 19931; iNordstrom et al. 


|20o4 


Reid et aZ.ll2007l; 


iHolmberg. Nordstrom & Andersen! 


2007) 


find that older 



stars tend to have lower metallicities on average. As it is 
difficult to estimate the ages of individual stars, the precise 
form of the relation is still disputed. However, there seems to 
be general agreement that there is a spread of metallicities at 
ea ch age, which is also su pported by other studies (Che n et 
aZ. l2003l;lHavwoodll2008l ; IStanghellini fc Havwoodll2010h . A 
metallicity spread amongst coeval stars is inconsistent with 
a simple chemical evolution model in which the metallicity 
of the disc rises monotonically in each annular bin, without 
mixin g in the radial direction. Again Sell wood fc Binnevl 
(2002) and Roskar et al. (2008ab) showed that when the 
disc supports recurrent transient spirals, the needed radial 
mixing arises naturally through angular momentum changes 
at cor otation that do not heat the disc. ISchonrich fc Binnevl 
(2009) developed the first chemical evolution model for the 
Milky Way disc to include this radial churning. 

The transient nature of the patterns is an essential as- 
pect to produce efficient churning. The horse-shoe orbits 
near co-rotation of a large-amplitude spiral cause a single 
change in the angular momentum if the spiral grows and 
decays in less than half the (long) orbit period in the rotat- 
ing /rameQ Were the spiral to persist for longer than this, 
then the exchanges a star experiences at one arm would be 
undone when the star encounters the next arm, and the star 
would merely take periodic, equal inw ard and outward steps 
in radius leading to no lasting change. iMinchev etai\ (|2010h 
report substantial mixing due to combined influence of a bar 
and spirals, but their simulations were of short duration and 
it is likely the effect is large only as the bar forms. Only re- 
current transient patterns, with corotation radii spread at 
random locations, can cause the home radii of stars to "dif- 
fuse" across the disc throughout its life. 



1 In fact. lSellwood fc Bimievl (2002) argue that the spiral mode 
grows until it reaches the amplitude at which the periods of horse- 
shoe orbits become short enough for this to happen for many 
stars, which causes the disturbance to disperse. 



6.4 Discussion 

Tremaine (private communication) pointed out that tran- 
sient spirals could co-exist with long-lived patterns. This 
possibility could be consistent, for example, with the change 
in the appearance of sp iral patterns between the blue and 
ne ar-IR passbands (e.g. iBlock fc WainscoatJ 199ll; Block et 
al. 1 1994t h although these data contain no information about 
the lifetimes of features of any colour. The co-existence of 
short-lived spirals seems inconsistent with th e specific mode 
theory for long-lived patterns proposed by iBertin fc Linl 
(1996), since transient waves would heat the outer disc to 
an extent that must render their proposed mechanism un- 
workable, as shown in However, some other possible (as 
yet unknown) mechanism for long-lived spirals might remain 
viable. 

An observational test for the coexistence of both short- 
and long-lived patterns would be a daunting challenge. For 
a well-constructed sample of galaxies, one would have to de- 
termine the mean amplitude needed to produce the desired 
heating and churning by the short-lived patterns, and then 
show that the mean observed spiral amplitudes, estimated 
somehow from the passband dependent photometry, would 
or would not allow an additional significant long-lived spiral 
component. 



7 CONCLUSIONS 

The question of the lifetimes of spiral patterns is impor- 
tant because galaxy discs evolve to a much greater extent 
if spirals are transient than if they are quasi-steady (jj2j. In 
this paper, I have summarized the evidence that bears on 
the question; were any one piece decisive, there would have 
been no need to write this paper. However, I find that the 
transient picture is favoured by the velocity structure in the 
solar neighbourhood ( 93J> and by the behaviour of simula- 
tions (§Sf4l4 &[5]). 

Direct observational evidence from external galaxies 
(32| is frustratingly inconclusive, largely because our sin- 
gle snapshot view can tell us little about the lifetimes of 
the patterns or, if self-excited, the mechanism that caused 
them. However, the overall appearance of the velocity distri- 
bution of stars in the solar neighbourhood is most naturally 
accounted for by a succession of transient spirals, with indi- 
vidual features perhaps being attributable to resonances. 

Theoretical work (iJJ) over many decades has pursued 
the two mutually exclusive views. In their efforts to gain 
the upper hand, the rivalry has spurred great improvements 
in our understanding of density wave mechanics. However, 
neither side has been able to land a knockout blow. 

Simulations have long manifested short-lived transient 
features, and I have presented two new results here. The 
first, in §U4, shows that galax y models with cool, low-mass 
discs of the type proposed bv lBertin et all (|l989l ) as likely 
to manifest long-lived, bi-symmetric spiral patterns will in- 
stead heat quickly as a result of short-lived transient spiral 
patterns with m > 2. While the behaviour in the simu- 
lations is not fully understood, swing-amplification gives a 
clear reason for the preference for higher sectoral harmonic 
disturbances in low-mass discs. Naturally, my demonstra- 
tion that the most carefully worked-out theory for long-lived 
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modes may not be viable, does not exclude the possibility 
of long-lived modes altogether. 

I have also shown (£01) that the most clearly presented 
case of a long-lived spiral in a simulation is not, in fact, a sin- 
gle spiral mode, but a sequence of a few strong patterns that 
each have short lifetimes. Thus I am unaware of a credible 
case of a long-lived spiral in any simulation, and therefore 
simulation results provide strong support for the transient 
hypothesis. The principal reason this is not decisive is that 
we still do not fully understand the origin of spiral patterns 
in the simulations. 

Indirect evidence (§6$ suggests that three quite diverse 
properties of galaxies are naturally explained by transient 
spiral waves, but an underlying long-lived wave is a possi- 
bility that would be very difficult to exclude observationally. 

To hold that spirals are long-lived, one must also argue, 
(1) that significant long-lived spirals underlie the transients 
that appear to be needed to explain the velocity structure 
in the solar neighbourhood, the high random motions of old 
disc stars and radial mixing in discs, (2) that the long-lived 
pattern adjusts its rotation rate slowly over time, so that 
fresh stars become the sources and sinks of angular momen- 
tum to drive the pattern, and (3) that no simulation has 
yet been able to capture the correct dynamical behaviour of 
spiral patterns. 

However, the entire question of spiral lifetimes will not 
be settled until a theory of spiral pattern generation, which 
is supported by observational data and consistent with sim- 
ulation results, gains wide acceptance. 
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